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We present a single band Hubbard model on honeycomb lattice with link-symmetric spin-orbital 
interaction. In the large U limit, this model is mapped to a conventional antiferromangentic (AFM) 
Heisenberg model combining with a 'zig-zag' AFM Heisenberg one. The latter AFM order which 
comes from the symmetric spin-orbital coupling was favored by a recent experiment for A2lr03 
(A=Na, Li) [6-9 . Competition between two different AFM orders leads to a spin liquid phase 
characterized by Kitaev model that happens if the spin-orbital coupling has the same strength as 
the hopping amplitude. This deepens our understanding of the origin of the Kitaev coupling: It 
comes from the balance between the hopping and spin-orbital coupling but to get rid of the hopping. 
In a moderate U, this model has a charge-spin gapped phase with two different spinon excitations. 
In a weak interaction, the original gapless Dirac fermion is gapped by the spin-orbital coupling while 
a new gapless Dirac fermion appears at the same Dirac points. 

PACS numbers: 71.20.Be, 75.25.Dk, 75.30.Et, 75.10.Jm 



Kitaev model on honeycomb lattice has attracted a 
lot of interests because of its exact solubility, Majo- 
rana fermion excitation, non-trivial topological orders 
and amazing non-abelian anyon pQ. This non-abelian 
anyon is the key object to design a device of topological 
quantum computer 2J. A possible way to realize Ki- 
taev model on optical lattice has been proposed [3]. Re- 
cently, studies of iridium oxides showed that this peculiar 
quantum spin liquid (SL) is possible to be realized [U [5] ■ 
While the experiments for A2lr03 (A=Na, Li) supported 
a 'zig-zag' antiferromagnetic (AFM) phase [SHS], Kitaev 
SL is possible to be observed in a high pressure [7] . These 
progresses drive a number of further researches |10H17j . 

It is argued that the strong spin-orbtial(SO) inter- 
action plays an important role in deriving the Kitaev- 
Heisenberg (KH) model |3]. However, the exact soluble 
Kitaev model is on the opposite side of the Neel order, 
i.e., one has to get rid of the hopping to obtain it. 

In this Letter, we will directly show how the SO in- 
teraction induces a Kitaev coupling in the context of a 
single band repulsive Hubbard model. What we observe 
is as follows: For a large U, while the hopping between 
the nearest neighbor sites contributes to the Neel order 
of spin S, the SO coupling offers a 'zig-zag' AFM order 
with a variant spin S. The latter one is the most likely 
spin structure in A2lr03 according to the recent experi- 
ments [6-9] . The competition between these two different 
AFM orders gives birth to a SL phase delegated by Ki- 
taev model with anyonic excitations. In fact, when the 
SO coupling strength is the same as the hopping ampli- 
tude, the resulting spin model is exactly Kitaev model. 
This is the crucial difference of our result from that in the 
HK model: The Kitaev phase is in the midway between 
the hopping amplitude and the SO coupling but getting 
rid of the hopping. We notice that the Kitaev phase may 



also appear in the competition between two different FM 
phases . 

The SO interaction which can lead to previously men- 
tioned physics is a time-reversal invariant SO interac- 
tion with the nearest neighbor site, i.e., link-dependence. 
The SO couplings with an orientation-dependence ex- 
ist in many Mott insulators with orbital orders and 
fluctuations [18 . In this Letter, we include a particu- 
lar link-dependent SO coupling to the Hubbard model 
and call it the symmetric SO-Hubbard (SSO-Hubbard) 
model. The strong coupling Mott physics has been sum- 
marized in the previous paragraph. For a moderate U, a 
gapped SL phase has been predicted even for the Hub- 
bard model without a SO coupling[T5]. A slave boson 
theory [201 HI] also showed that there is such a SL state 
at the mean field level but may be broken down due to 
a staggered U(l) gauge fluctuation for some symmetries 
are broken [52]. This SL is also found in a Kane-Mele- 
Hubbard model [53]. What we find is that, by using the 
slave boson method [32], the SSO coupling will enlarge 
the original spinon gap while a new branch of spinon ex- 
citation takes the gap of the original one. In the weak 
coupling limit, we also find that the well-known Dirac 
fermion on the honeycomb lattice is gapped by the SSO- 
coupling while a new gapless Dirac fermion whose veloc- 
ity is solely determined by SSO coupling is exactly locates 
at the original Dirac points. 

We consider the SSO-Hubbard model on honeycomb 
lattice with links a = x,y,z described by the following 
Hamiltonian 

Hssoh = H T + H sso + H v = -t^2{c\ a c ]a + h.c) 

a;(ij) a 1 



2 



where Vij — —Vji = 1. Ci a is the electron annihilation 
operator at site i. a a are Pauli matrices. U > is 
Hubbard potential, (ij) labels the summation over the 
nearest neighbor sties and (ij) a means over all the a- 
links. Different from that in Kane-Mele model [2"1] , here 
the spin orbital coupling is x, y, z-cycling symmetric and 
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where J = a ^= b ^ c and S = ^c^Uca'Ca 1 is the elec- 
tron spin operator. At the half filling, the SSO-Hubbard 
model reduces to 

Hjj = Jj^ Si • Sj + J ( S i S ? - S i S j - S Z S t)> ( 3 ) 

(ij) a;{ij) a 

where J = as well-known. The crucial difference here 
from that in the KH model is both J and J are positive, 
which means we have two different AFM orders. For 
A = 0, it is AFM Heisenberg model. If t = 0, it is another 
quantum liquid which differs a sign from the stripy AFM 
model and is of a 'zig-zag ' AFM order (See Fig. [TJ. It 
seems to be the most likely spin order in A 2 Ir03 in recent 
experimental reports combining with the first principles 
calculation [SHS] or a direct measurement [5]- This 'zig- 
zag' AFM model in fact can be mapped to an AFM model 
Hj = Jj2(ij) Sj • Sj with Sj defined as follows: Divide 
the honeycomb lattice into four sublattices, keep Sj = Sj 
in one sublattice and change in other three, e.g., Sj = 
(Sf,—Sf,—Sf) for the cc-link direction sublattice (and 
cycling for the y- and z-link's. See Fig. fija) [3].) 




(a) Stripy AFM (b) Zig-zag AFM 



FIG. 1: (Color online) (a) The stripy AFM order. Four 
different colored spots label four sublattices. (b) The zig-zag 
AFM order. The black and white spots divide the lattice in 
to a bipartite one. 

Rewrite Eq. ([3| as 

Hjj = {J- J) E S < ' S J + 2j E S i S j ( 4 ) 
{ij) a-{ij) a 

= {J-J)^S i -S j + 2J E SfSf- (5) 

(ij) a;{ij) a 



involves in the nearest neighbor sites instead of the next 
nearest neighbors. 

Strong coupling limit. We consider Ht and Hsso as the 
perturbations. Besides the well-known AFM coupling as- 
sociated with Ht, the Hsso contribution, to the second 
order for a given link a , is given by, 



I 

One immediately sees that if t — A, we arrives at Kitaev 
model with Jk = 2 J. It is coaction both of the hopping 
and the SSO coupling gives birth to the Kitaev coupling. 
Notice that the system, either for J = or J = 0, is in the 
conventional Landau symmetry breaking ordered phase. 
The Kitaev phase, however, is a SL phase without any 
Landau-type order parameter. That is to say, the Kitaev 
physics is the result of competition between two different 
AFM orders. There are quantum phase transitions from 
the Neel phase to the Kitaev SL and then to the 'zig- 
zag' AFM phase. The ground state energy per site for 
the Neel order is estimated as ejy/J w (1 — /3/3)(Sj • 
Sj) for f3 = J I J and (Sj • S i ) -0.37 0] while the 
upper boundary of that of the Kitaev SL is e^/J < (1 + 
P){S?Sf) for (SfSf) = -0.13 [25]. This estimates the 
critical point /3 C < 0.95. According Eq. another 
critical point is l/j3 c « 1.05. By exact diagonalization 
calculations, up to 24 lattice sites, we confirm the critical 
point is very close to the above estimate but are unable to 
determine that according to the ground state properties. 
Since the gapless excitation in the Kitaev phase is very 
different from that in the AFM phases: The former is 
Majorana fermion while the latter is bosonic spin wave. 
Therefore, we can estimate the critical points by looking 
for the level crossing of the first excitation state energy 
levels of two different phases in a finite size (See Fig. H 
(a)). (We also use this criteria to check the critical points 
in the KH model which are the same as the results in 
Ref. [IJ.) The 24-site data gives f3 c sa 0.93 for the phase 
transition from the Neel order to the Kitaev phase. 

In Kitaev model, if a perturbed external magnetic field 
is applied, the SL becomes gapped and the phase we are 
discussed turns to be topologically non-trivial with a non- 
vanishing Chern number and non-abelian anyons. We 
expect these non-trivial properties exist in the whole SL 
phase and a gapless chiral excitation emerges in the edge 
of the system. Both of the AFM states are also gapped. 
At the critical points, the gaps collapse. 

When the system deviates from the half filling, one can 
study the doped t-X-J-J model. 

Original Kitaev model allows a phase transition from 
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FIG. 2: (Color online) (a) The level crossing of the energy 
levels of the first excited states for lattice sites N = 8, 16 and 
24. (b) Schematic phase diagram in J-J plane. 



the non-abelian phase to three equal gapped SL phase 
with abelian anyon. This can also be realized in this SSO- 
Hubbard model by introducing link-dependent hopping 
amplitudes t x ,t y ,t z and SO couplings X x ,X y ,X z . In this 
way, we can have a whole parameter Kitaev phase as 
well as the fruitful phase diagram from the XYZ model 
to Kitaev model and that in Kitaev model. 

In the KH model studied in Ref. [4], there is a phase 
transition from the stripy AFM phase to the Kitaev 
phase. This can not be explained by the single band 



SSO repulsive Hubbard model. Instead, we find that it 
happens when two different FM orders |26| compete, i.e., 
if J < and J < in Eq. one has 



H 'jj 



(\J\ 



Qa Qa 



(6) 



(Or expressed by S with J O J as Eq. |5|.) When J = 0, 
one has the FM order while it is the stripy AFM order [J| 
for J = 0. Again, Kitaev model comes out at J = J. The 
phase transition from the Kitaev phase to the stripy AFM 
phase happens at J = 2 J/3 according to the result by 
Ref. [1]. Thus, the window of the Kitaev phase is larger 
than the AFM case: 2/3 < J/J < 3/2. This captures 
the part of the KH model [3] from the stripy AFM to the 
Kitaev phases. Furthermore, the phase transition from 
the AFM Neel order to the stripy AFM order |4j is the 
result of the competition between J > and J < 0. A 
model with J < and J > corresponds to a phase 
transition from the FM to the zig-zag AFM order. A 
schematic phase diagram is depicted in Fig. [2jb) . 

Weak interaction limit. When U = 0, the free-electron 
gas on honeycomb lattice has a gapless Dirac fermonic 
low energy spectrum in two valleys near the Dirac points. 
Now we discuss the effect of the SSO coupling to the 
Dirac semi-metal. The honeycomb lattice is bipartite 
to A and B sublattices (See Fig. [TJb)) . Using the z- 
link as the elemental cell, the electron annihilation op- 
erator may be labeled by the sublattice suffices, i.e., 
Cut — > {cAa-,i, CBa,i)- Under the Fourier transformation 
c aa = ^ k c QCT ke 4k ri , the Hamiltonian with the SSO 
coupling is given by 







H T + HsSO = E( C Atk> c Bfk> C W> C i 
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where f(k) = 1 + e lk ni + e ik n2 and g±(k) = i(e lk ni ± 
ie ik n2 ) with the vectors rix = (1/2, \/3/2) and n 2 = 
( — 1/2, \/3/2). The spectra of quasi-electron excitations 
then can be easily calculated: 

E\ = (6A 2 + 2< 2 |/| 2 )/2 (8) 

±V(6A 2 + 2*2|/|2)2 - 4|A2 3+ . g _ - i2/ 2 - A 2 | 2 /2 

When A = 0, we recove E = ±t\f(k)\ cx q near Dirac 
points K± = (±4tt/3, 0) ( q = k — K±). Now, because 
|/| 2 oc q 2 and \g+(k)g-(k) — 1| cx q near K±, we have 

E- k ±X\g+g- - 1| ~ ±Xq, E+ = ±V&X + 0{q 2 ). (9) 



The i-dependent branch (E+) of fermnionic excitation is 
gapped while a new gapless Dirac fermion mode (£L) 
with a speed A appears at the same Dirac points. This 
Dirac fermion differs from the free one as the latter does 
not mix the spins and the former mixes the spins. 

Moderate interaction. In a moderate U, it is predicted 
that there is a featureless gapped SL phase [19]. A 
slave boson study showed that there is also a charge- 
spin gapped phase in a moderate U region [55] although 
the magnitude of this region is numerically different that 
from the quantum Monte Carlo study in Ref. [TH] . More- 
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over, the SL state is broken down due to translation 
symmetry and lattice rotation symmetry breaking from a 
staggered U(l) fluctuation. As we have seen before new 
physical properties of the system arise from the SSO in- 
teraction. We now use the salve boson technique to see 
what is the effect of this SSO coupling to the system in 
this interaction region. 

We briefly review the slave boson theory (20] |2T] , 
There are four different states in a single site: |0),| t 
), | |), | U) with the completeness |0)(0| + | t)(t I + 
I + 1 ti) (ti 1 = 1- The electron operator can be expressed 
as c a = |0)(er| +er| — cr)(t| |. At half- filling, one can thing 
each site has averagely an electron charge. If subtract- 
ing this average charge, the state |0) can be thought as 
carrying a positive charge while | fi) carries a negative 
charge and \a) carries spin degree with charge neutral. 
In the operator formalism, one has c a = f a + aj^_ a d. 
This is the slave boson technique with a spin-charge sep- 
aration if ~ |0),d^ ~ | t-i) i called holon and doublon, 
are bosonic operators carrying charge ± while /J ~ \a) 
is called charge neutral fermonic spinon. The complete- 
ness condition becomes a constraint on these charge and 
spin degrees: S* = n hl + n di + J2 a n f^i - 1 = for every 
site. Within the slave boson theory, then, the effective 
Hamiltonian is given by 

Hefj = uj24di-Yl E [(**/,« - x xhj)x b ,ij 



The eigen energies then are easily obtained by 



a a— links 



(tA£y - AA%.)A fcjl + h _ c] +iJ2^iSi, (10) 



where x%j = E CT flfjv, Xf,ij = ^'^ijfL<a'fj<r, 
X6,« = h\h 3 - d\d r And A^. = £ o rr f, t// , t . A a f>ij = 

Y.<ja> a ' v fj fi° (J aa' fj-<y' and ^b,ij = h t dj +dihj. For a 
non-zero A/ then A&, the pair nature of the boson wave 
function and the holon-doublon symmetry require \b = 
and then tx° f - \x a f = 0. 

Now we follow the same trick used in Ref. [55] to solve 
the mean field theory. The boson problem in our model 
is not much different from that in Ref. [22 if tA° in 
Ref. [55], the mean field value of A° i:) , is replaced by 
t A" — AA^ in every link. The charge excitation spectrum 

is E± = ±| + ^(^0)2 _| A ^ where £ = (if) and 

A fk = J2 a ( tA °f ~ AApe lk n " for n x ' y = m i2 and n z = 
[55] , The mean field effective Hamiltonian for spinon is 
given by 



H f = -(V2)^(/L,jA^ + /L, fe w) (ii) 

k 

k 

k 

fkkfk-k X 9 + (k)]A b + h.c. 



E 2 = Co/ 4 +[6A 2 + 2P |/| 2 



(12) 



±(1/2W(6A 2 + 2^|/|2)2 _ 4 |A 2 5+ ^ - i 2 / 2 - A 2 |^ 



where t — tAh 



and A = A Af,. 

■2 , 



Notice that in Eq. (12 1, 



except the first term £q/4, the rest is of the same form 
as Eq. ([8]), that of the Dirac fermions at U = 0. There- 
fore, there are two branches of spinon excitations whose 
energy minima are at the Dirac points and the gaps are 

|^o 1/2 and y^/4~T 6A 2 , respectively. As the same as 
that in Ref. [35] , our mean field states also suffer a stag- 
gered U(l) global gauge fluctuation as fiA — > e la fiA and 
fiB —> e ~ la fiB and the SL state is possibly broken down 
in a semiclassical limit. Besides the next nearest neigh- 
bor hopping may destroy this U(l) gauge symmetry [22], 
the Kane-Mele SO coupling [21] also concerns the next 
nearest neighbor sites and breaks this U(l) to Z 2 gauge 
symmetry so that the system is turned to be a SL. More- 
over, while the Kane-Mele SO coupling gives a mass to 
the new Dirac fermion at the weak coupling and the sys- 
tem is expected to be a topological insulator with the 
time reversal svmmetrv|24j. the strongly-corrected be- 
haviors with this SO coupling are more intriguing |27j . 
Finally, if the system is doped, those SL phases may be 
superconductors and the competition between the d-wave 
superconductivity and a p-wave superconductivity is ex- 
pected as U varys. 

In conclusions, we have studied the physical behaviors 
of the SSO-Hubbard model in the strong, weak and mod- 
erate repulsive interactions. We showed that a SL phase 
can be a result of a competition between two Landau- 
type ordered states: the Kitaev phase appears when two 
different AFM couplings become comparable, i.e., both 
the hopping amplitude and the SO coupling play equal 
roles. Meanwhile, the hopping does not play any role in 
the limit of U — > 0. A phase with both charge and spin 
gapped emerges in the midway between the weak and 
strong couplings. We do not determine the critical U for 
these different phases. Since our mean field equation for 
the slave boson's is formally the same as that given in 
Ref. [52J, the result is the same. The critical interaction 
into the AFM orders is also very close that when A = 0. 
We have seen that the Kitaev phase appears due to the 
competition between two different AFM orders. This im- 
plies the Kitaev phase must not appear before these two 
AFM phases but once these two ordered ones exist, there 
should be this SL as a result of the competition. 

In conclusions, we have studied the SSO-Hubbard 
model for the repulsive Hubbard U. This work was sup- 
ported by NNSF of China and the 973 program of MOST 
of China. 
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